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We investigate the figure of merit of a quantum dot (QD) in the Coulomb blockade regime. It is found that the 
figure of merit ZT may be quite high if only single energy level in the QD is considered. On the other hand, with 
two or multi energy levels in the QD and without the Coulomb interaction, the ZT is strongly suppressed by 
the bipolar effect due to small level spacing. However, in the presence of the Coulomb interaction, the effective 
level spacing is enlarged and the bipolar effect is weakened, resulting in ZT to be considerably high. Thus, it is 
more likely to find a high efficient thermoelectric QDs with large Coulomb interaction. By using the parameters 
for a typical QD, the ZT can reach over 5. 
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I. INTRODUCTION 

Thermoelectric materials are such materials that can di- 
rectly convert the thermal energy into the electrical energy. 
Thus the thermoelectric energy conversion technology has 
been recognized as the most feasible energy conversion tech- 
nology. However, due to its low efficiency, this technology 
has not been widely used. Thus it is important to find high ef- 
ficient thermoelectric materials. The efficiency of thermoelec- 
trical materials is measured by the dimensionless thermoelec- 
trical figure of merit ZT, while ZT = gS 2 T/h. Here S is the 
Seebeck coefficient, a is the electric conductivity and k is the 
total thermal conductivity which contains the lattice thermal 
conductivity k; and electric (carrier) thermal conductivity K e , 
and T is the operating temperature of the device. 1 For a mate- 
rial to be a good thermoelectric material, it must have a large 
ZT, which means in order to achieve a large ZT, one must 
increase the Seebeck coefficient S and electric conductivity a 
and decrease the thermal conductivity k. However, it seems 
difficult to have a high ZT in nature materials. Several rea- 
sons hinder the rising of ZT. First, in conventional solids the 
Wiedemann-Franz law (K e /aT = (fcs7r) 2 /3e 2 ) is obeyed, 2 
which means that an increase in the electric conductivity also 
leads to an increase in the thermal conductivity. Second, ac- 
cording to the Mott relation, 3 an increase in the electric con- 
ductivity is apt to lead to an decrease in Seebeck coefficient. 
Thus in the past fifty years, the maximum ZT is holding at 
about 1, This largely affect the industrial applications. 

Recently, the advance in nanostructure materials have 
largely stimulated the development in thermoelectric materi- 
als. Due to the quantum phenomena emerged in nanostructure 
materials, the classical results such as the Mott relation and 
the Wiedemann-Franz law may not hold. 4 What's more, the 
thermoelectric properties of the nanostructure materials can 
be modulated by changing the gate voltage. Thus it opens 
a new and wide road to find efficient thermoelectric mate- 
rials. The idea of using low dimensional structure materi- 
als to gain high ZT was first introduced by Hicks and Dres- 
selhaus in 1993. 5 They theoretically show that ZT increases 
swiftly as the dimensions decrease, far beyond the value ob- 



tainable in bulk materials. Following this suggestion and with 
the development in nanotechnology, various groups were able 
to fabricate nanostuctures and measure their thermoelectric 
properties. 6-10 For example, Harman et. al. have measured 
the thermoelectric properties of quantum dot and a maximum 
value ZT w 2 was obtained. 7 Venkatasubramanian et. al. 
have measured a thin-film thermoelectric device and have ob- 
served a maximum ZT of ~ 2.4 at room-temperature. 8 Apart 
from the experimental efforts, many theoretical studies have 
been carried out on low dimensional structures such as quan- 
tum dots, nanowires, and superlattices^i^i 4 . For example, 
Venkatasubramanian and Chen have concluded that the main 
reason of high ZT in low dimensional materials is due to a 
significant reduction in lattice thermal conductivity. 13 A giant 
figure of merit in single-molecule device is obtained by Finch 
et. al. J^- All these efforts show that a high figure of merit may 
exist in nanomaterials. However, due to the complexity and 
expensiveness of the nanomaterials, there still a long way to 
go for the commercial applications of nanostructure thermo- 
electric materials. At present the most promising nanostruc- 
ture thermoelectric material is nanocomposite thermoelectric 
materia h 2 ' 15 ' 16 

In this paper, we study the thermoelectric properties of 
lead-QD-lead system with the QD in the Coulomb block- 
ade regime. The thermoelectric properties of QD have been 
widely studied. For example, Beenakker et al. have inves- 
tigated the thermal properties of QD with multiple energy 
levels, 417 but they only considered the situation that temper- 
ature (fcflT) is much bigger than the level-width (T). On the 
other hand, Turek et al. and Murphy et al. have studied the 
thermal properties of the QD in the situation fc^T ~ r^S^ 
but the QD considered contains only a single energy level. 
Here we consider the QD having the multiple energy levels 
and in the situation of the temperature smaller than the in- 
teraction U. In this regime, some new phenomena emerge. 
By using the Landauer-Biitticker formalism combining the 
nonequilibrium Green's functions, 20 the electronic conduc- 
tivity, Seebeck coefficient, and thermal conductivity are ob- 
tained. Due to the electron-hole symmetry, the Seebeck coef- 
ficient is always antisymmetric. If only a single energy level 
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in the QD is considered, the ZT increases monotonously with 
temperature T, and ZT can be very large at high temperature, 
consistent with previous results. However, with two or multi 
levels and the temperature is on the order of energy gap, the 
ZT is strongly suppressed by the bipolar effect, mainly caused 
by the antisymmetric property of the Seebeck coefficient. On 
the other hand, when the Coulomb interaction U is consid- 
ered, the energy spacings are enlarged due to the Coulomb 
blockade effect. The bipolar effect is greatly reduced and 
high value of ZT may again be achieved. In a typical QD, 
the Coulomb interaction U is usually larger by an order than 
the linewidth T and the single particle energy spacing. Under 
these conditions, the ZT can be quite high, with its maximum 
value reaching over 5. 

H. MODEL AND FORMALISM 

The system of the lead-QD-lead can be described by the 
following Hamiltonian: 

H = Sakcl k(J C a ko- + t a k{d\ a C a k<j + H.C.) 

x U x 

+ 2^ £iTiia + ~2 TlicrTljcr' i (1) 

£=1,2;<7 i,a t j,a' {ia^jcr' ) 

where n irj — d\ a d ia , and a — L,R represent the left and 
right leads. c' aka and d\ a create an electron with spin a in the 
a lead and the ith energy level of QD, respectively. Here the 
intra-dot electron-electron Coulomb interaction is considered, 
with the interaction strength U. The second term in Eq. ([TJ 
describes the tunneling coupling between the QD and the two 
leads and t a k is the hopping matrix element. 

By using nonequilibrium Green's function methods, the 
electronic current and electric thermal current from the left 
lead flowing into the QD can be written in the forms: 20 ' 21 

where f a — f(oj — /i Q ) = l/{exp[(u) — /i q )//cbT] + 1} is the 
Fermi distribution of the a lead and T(w) is the transmission 
coefficient. T(lu) can be expressed by the following expres- 
sion: 

TH=Tr[iiJf (G r -G«)]. (3) 

Here G r (G a ) is the standard retarded (advanced) Greens 
function of the QDr ^ and T a ^ = J2 k ^\t a k\ 2 5{uj - e ak ) 
is the linewidth functions which assume to be independent 
of the energy uj. We introduce the following integrals I n (T) 
(n = 0, 1, 2, ..): J„(T) = -(2/h) J u n (df /8w)T(w)dw. By 
using the quantities I n (T), the linear electric conductance G, 
thermopower S, and thermal conductance k can be expressed 
in very simple forms:— 

G = e 2 / (T), (4-1) 
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FIG. 1: (Color online) G{2e 2 /h), k(2k B /h), S{k B /e), ZT vs. the 
level ei for the different temperature T for the single-level QD and 
(7 = 0. 

S= -h(T)/[eTI (T)], (4-2) 

k = (1/T)[/ 2 (T) - lf(T)/I a (T)}. (4-3) 

Therefore, the only question left is to calculate the Green's 
functions of the QD. 



III. NUMERICAL RESULT 

In the numerical investigation, we consider the symmet- 
ric barriers with Fl^j = ^R,ij = I\ and set T = 0.5 and 
(in = as the energy zero point. We consider the linear 
regime, then /i^ = fi^ = /i = 0. First, we study the 
case with QD possessing only one energy level £\ and in 
the absence of the Coulomb interaction (U = 0). In this 
case, the Green function of QD can be easily obtained as 
G lcr (w) = Gf*(w) = l/(u; — ei+«r). Inserting this into Eq. 
(0, the transmission coefficient T(uj) can be obtained and the 
thermoelectric properties can then be calculated straightfor- 
wardly. Fig. 1 shows the electric conductance G, the thermal 
conductance k, the thermopower S, and ZT versus the level 
E\ for the different temperature T. Variation of £i is equiva- 
lent to variation of the gate voltage in an experimental setting. 
The electric conductivity G and the thermal conductance k 
exhibit a single resonant peak at the position E\ = 0. The 
peak height of the thermal conductance k first increases and 
then decreases while the increase of the temperature ksT (see 
Fig. 1(b)). The reason is as follows. The thermal conductance 
is determined by two aspects: the heat transferred by each 
electron and the tunneling probability of each electron. When 
temperature increases, the average tunneling probability de- 
creases but the heat transferred by each electron increases, 
leading a non-monotonic relation of k and ksT . The prop- 
erty of thermopower S is described in Fig. 1(c). Here we can 
see that the curves are antisymmetric due to the electron-hole 
symmetry. The reason is as follows. The thermoelectric ef- 
fect is caused by the temperature difference. There are more 



FIG. 2: (Color online) G(2e 2 /h), k(2k B /h), S{k B /e), ZT vs. V g 
for the different temperature T for the two-levels QD and with the 
parameters U = and level interval Ae — 4. 



FIG. 3: (Color online) G(2e 2 /h), k(2k B /h), S(k B /e), ZT vs. 
V g for the different temperature T for the two-levels QD with the 
parameters U — 50 and level interval Ae = 4. 



electrons being excited above the chemical potential fi in the 
hotter region and correspondingly more holes being gener- 
ated below fi. When the energy level of QD is below fi, the 
main carriers are holes and then the thermal power is posi- 
tive. When the energy level is above /i, the main carriers are 
electrons and thus the thermal power is negative. So one can 
adjust the gate voltage or equivalently E\ and obtain the op- 
timized thermal power. Once the thermal power, the electron 
conductivity and the thermal conductivity are known, ZT can 
be calculated. Fig. 1(d) describes ZT as a function of QD's 
level £\ . The optimized ZT can be obtained by modulating E\ 
when the system is kept at a fixed temperature. With increase 
of temperature, the value of optimized ZT also increases and 
it goes to infinity as T approaches infinity. This is consistent 
with the previous result. 24 Of course, this is a non-sensible re- 
sult due to only single level being considered here. As we can 
see, at temperature T = 1 (i.e. T = 2T), the optimized ZT is 
about 2.5. 

In a realistic situation, depending on temperature, dot size, 
etc., one normally has to consider multi levels. For a multi- 
level dot, the spacing between neighboring levels is an impor- 
tant quantity. For that purpose, investigating a two-level dot 
will capture essential physics due to the level spacing. In the 
following we study the thermoelectric properties of a QD con- 
taining two energy levels. The Green function of QD with two 
energy levels is: 
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Here, the boldface letters (G, g, and X) represent the 2x2 
matrix. gJJ. is the Green Function of QD without coupling to 
the leads, g£. can be obtained from the equation of motion 
technique (the detailed deduction can be seen in Appendix A) 
and S r < are obtained from Dyson equations ( Here we just 



consider the first order of self-energy correction and have ne- 
glected the higher order of self-energy correction that due to 
the e-e interaction )jS 
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where N a = rii S + n- ia + n- is , n J(T is the electron occupation 
number in the ith energy level with the spin state a, a =\. 
while a =j" and a =t while a =\., and i = 1 while i = 2 and 
i = 2 while i = 1. In Eq.(|7J, [Nia] means the integer part of 
Ni a , {No-} = No- — [N a ], namely the decimal part of N a - 
In addition, the electron occupation numbers ni a need to be 
self-consistently calculated with the self-consistent equation 
ria = -if(du/2ir)G<(u). 

Fig. 2 shows the conductance G, the thermal conductivity 
k, the thermopower S, and ZT versus the gate voltage V g in 
the absence of the Coulomb interaction (U = 0). Here the 
energy levels and the gate voltage are related by: e\ = V g 
(5) and £2 = V g + Ae, where Ae is the spacing between the two 
levels. At the low temperature, G and k in Fig.2(a,b) exhibit 
two peaks due to the two energy levels. There seems no great 
change in G in comparison with that for the single-level QD. 
However, behavior of k in Fig. 2(b) is more sensitive to tem- 
perature. When temperature is of the order of Ae, the peaks 
are broadened to a degree to give rise to a huge peak. Mean- 
while, the thermopower S and the ZT are largely suppressed 
when the lead's chemical potential [i is between the two en- 
ergy levels (£2 > /i > ei). On the other hand, the optimized 
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ZT remains considerably large when fi is outside the two en- 
ergy levels (/i < E\, £2 or fi > e 1; e 2 ). In a real system there 
are many energy levels in a QD. The ZT with the chemical po- 
tential outside of the two levels is influenced by other levels. 
So we only focus on the ZT for the case with £2 > fi > £1, 
in which the optimized ZT is rather low. This is because the 
electron and holes are excited in the range of ksT, then the 
carriers can tunnel through the ith energy level in QD when 
fi — ksT < E.- L < fi + kgT. When temperature is significantly 
lower than the level spacing As, the carriers can only choose 
one energy level to tunnel, which is similar to the single-level 
QD. When temperature is on the order of Ae, the carriers can 
choose both levels to tunnel. Due to the temperature differ- 
ence between the two leads, the electrons above (below) the 
chemical potential fi in the left lead are more (less) than the 
electrons in the right lead. While e\ < fi < £2, electrons 
tunnel from the left lead to the right lead through the level 
£2, at the same time electrons tunnel from the right lead to 
the left lead through the level e\, in other words the hole tun- 
nels from the left lead to the right lead through E\. This is a 
bipolar effect: a nonzero heat conduction emerges even when 
the net electrical current is zero. Furthermore, the Seebeck 
coefficient S is significantly suppressed since the carriers are 
tunneling through the QD in both channels via opposite direc- 
tions. The above reasons cause ZT to be very small. 

From Fig. 2(d), we can see that the optimized ZT is arising 
with temperature when temperature is much lower than the 
level spacing Ae. However, when T is of the order of Ae, the 
bipolar effect is enhanced and ZT is decreased. In order to 
increase ZT, one needs to find a QD with large enough level 
spacing Ae such that to reduce the bipolar effect. In a typical 
QD, the bare value of Ae is not big enough. An alternative 
way is to choose a QD with strong Coulomb interaction U. 
Though the Coulomb interaction may suppress ZT^ at the 
same time it can broaden the energy level. The effective level 
spacings can be widened to U + Ae, which can give rise to a 
large ZT. In a typical QD, U can be quite large. For example, 
U in Cqq is on the order of 0.1 ev. 25 In addition, the region 
that reduces the bipolar effect is U > fcsT/10. In this situ- 
ation, the suppression of ZT due to Coulomb interaction can 
be neglected. 

Now we consider that the QD contains two energy levels 
and both inter-level and intra-level Coulomb interactions ex- 
ist. For convenience, we set the inter-level interaction to be 
equal to intra-level interaction. The interaction U is about one 
order larger than the level spacing Ae, as for a typical QD. 
In the existence of the interaction U, the levels are located at 
£1, £1 + U, £1 + Ae + 2U, and £1 + As + 3*7 due to the 
Coulomb blockade effect. Thus the linear electric conductiv- 
ity G exhibits four main resonant peaks at the positions of 
V g = 0, -U, -2*7 - As and -3*7 - As (see Fig.3a). In ad- 
dition, there are also some smaller peaks at low temperature 
due to the tunneling through the excite states. Fig. 3(b) shows 
the thermal conductivity n versus the gate voltage. We can see 
that at high temperature (ksT ~ U) the bipolar effect occurs 
and huge peaks of the thermal conductivity emerge at the val- 
leys between the two main adjacent peaks of the conductivity. 

The thermopower S shown in Fig. 3(c) is sensitive to the 




FIG. 4: (Color online) ZT vs. temperature T for the different 
coulomb interaction U with Ae = 4. For each situation we choose a 
suitable gate voltage. 

slope of conductivity. It is clearly seen that at low temperature 
the thermal power changes from positive to negative when the 
gate voltage moves across each peak of the conductivity (e.g. 
T=0.1). With temperature rising, the small peaks in the con- 
ductivity are absorbed to the main peaks and accordingly the 
transitions from positive to negative of the thermal power are 
reduced to four. The peak and valley values of the thermal 
power are enhanced with rising temperature. However, due 
to the bipolar effect, the peaks and valleys reach their maxi- 
mum values at about T = 5 (i.e. U/10). Further increasing 
of T will decrease these values. After G, n, and S are calcu- 
lated, the ZT can be determined, shown in Fig. 3(d). Since the 
Coulomb interaction U broadens the level spacing, the bipo- 
lar effect is greatly suppressed and ZT is enhanced. The opti- 
mized ZT can be over 5, much larger than the value in Fig. 2 
without interaction. 

To further investigate the effect of temperature and 
Coulomb interaction, we numerically calculate the ZT versus 
temperature at different Coulomb interaction U (see Fig. 4). 
Notice that the bipolar effect can be enhanced with increas- 
ing temperature but weakened with broadening of the energy 
spacing. Thus when fc^T < U/10, ZT is enhanced with 
rising temperature because of the weak bipolar effect. When 
ksT > U/10, ZT saturates with further increasing of tem- 
perature. Moreover, with increase of U, the optimized ZT 
also increases. 



IV. CONCLUSION 

We investigate the thermoelectric properties of a QD that 
contains one or two levels and is in the Coulomb blockade 
regime. The results exhibit that in the absence of the Coulomb 
interaction, the ZT can be very high if only one level in the 
QD is considered, but the ZT is greatly suppressed with multi 
levels due to the bipolar effect. When the Coulomb interaction 
U is considered in the QD, the spacings of energy levels are 
increased, and the bipolar effect is weakened, thus the ZT can 
be considerably high. For an actual QD in which its Coulomb 
interaction is one order larger than the level spacing, the op- 
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timized ZT can be over 5, much larger than the values from 
natural materials. 
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Appendix A: The Green Function of QD 



This is an exact solution without any approximation. How- 
ever, it is difficult to self-consistently calculate (ni a ni'cr') 
and (riiaTii' ernes) through numerical means, and some ap- 
proximations are needed. Here, we make the approximation 
(n ia n V(T i) = while (n ia ) + (rv CT /) < 1 and (n^rVo-') = 
(nia) + (nj/cr') — 1 while (rii a ) + (rv<r') > 1. In addition, 
we make another approximation (ni a riii a riii s ) = while 
N la < 2 and (n lS ni, a n. v 9 ) = {N irT } while N icr > 2. These 
approximations are reasonable since the fluctuation of the oc- 
cupation number in the QD is less than one at zero bias and 
the temperature fc^T < U. 23 Thus, when N ia < 1, the Green 
function can be simplified as: 

** = + u - ei -U }6 ^ (A3) 



In this appendix we give a detailed deduction on how to get when 1<N < % the Green function can be simp i ified as: 
the Green function of the QD with two energy levels. 

First, a single QD without coupling to the leads can be de- 
scribed by the following Hamiltonian: 



Hr 



U 



^ eih la + — ^2 h la h Jcr , (Al) 



i=1,2;ct 



i,c,j,<7-'(«7'^icr') 



For an isolated QD, the exact Green function can be obtained 
by the equation of motion technique: 



(u - Ej -U) + U((n lS ) + (r^y) + (7y a )) 

(u - £i)(uj -Si-U) 

2U 2 ((n l9 n l , cr ) + (n Va ni, s ) + (n lS n Va )) 
(uj - e t )(u - £j - U)(u - Si - 2U) 

QU 3 (n is n ilcr n lls ) 

(ui - ei)(u -ei- U)(u> -ei- 2U){uj -ei- 3U) 



= { 1-{AU + W } , 

x u-e, - U lo -e, - 2U { J ' 



(A4) 



When 2 < N irT < 3, the Green function can be simplified as: 



f l-{N ia } {N^} 

g - ^ { u-e l -2U + u - £i -3U }Sii ' (A5) 
Then the final form of Green function can be written as: 



;■» = { 



1 - {N la } 



{N la } 



u-Ei- [N ia ]U ' u - Ei - ([N i(T ] + l)^ 5 * 3 ' 

(A6) 
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